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Magnetic response of a superconductor depends on the thermodynamic stability of vortex in the 
material. Here we show that the vortex stability has a close relation with the ratio of the magnetic 
field at the vortex core center to the thermodynamic critical field. This finding provides a new 
criterion for the classification of the superconductors according to their magnetic responses. 
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Type-1 and type-2 superconductors exhibit different 
magnetic responses to externally applied magnetic field. 
Whether there exists stable vortex is the distinction be¬ 
tween two types of superconductors^. Due to the for¬ 
mation of stable vortices, mix state appears in a type-2 
superconductor under a certain range of external field. 
As a consequence of thermodynamic instability of vor¬ 
tex, only macroscopic Meissner state and fully normal 
state or the state of the coexistence of these two states 
(intermediate state) can exist in a type-1 material. 

To date, the identifications of types of superconduc¬ 
tors are based on the following three considerations: (1) 
to calculate the surface energy of a superconductor under 
the thermodynamic critical magnetic field^; (2) to com¬ 
pare the lower (higher) critical field at which the vortex 
entry into (exit from) the superconductor to the ther¬ 
modynamic critical magnetic field^^; (3) to analyze the 
interaction between vortices^—. All these three consid¬ 
erations show that, the magnetic response of a super¬ 
conductor is determined by a dimensionless parameter 
k, which is defined as the ratio of the penetration depth 
to the coherence length. And the critical value of the 
Ginzburg-Landau (GL) parameter k c = 1 /\[2 represents 
a boundary between two types of superconductors. 


tivity. Here we use the GL free energy^ 
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where /„o is the free energy density of the body in the 
normal state in the absence of the external field, e* and 
to* are the effective mass and charge of the Cooper pair, 
|\I/| is the modulus of the superconducting order parame¬ 
ter, B is the magnetic field, W(|| 2 ) = a |\I/| 2 + (3/2 |T| 4 . 
The free energy CD is equivalent to the usual GL model 
in which the gauge potential A and the order parame¬ 
ter L are the functions to describe the superconductiv¬ 
ity, and B = V x A. We consider an isolate vortex 
in an infinite sample, and use instead of the variable 
r, the function |L| and B the dimensionless quantities 
p = r/A, \ip\ = I'Ll/W 0 , B = |B| /H c , where A = 
( 'm*c 2 (3/A'Ke* 2 |a|) is the penetration depth, Lo = 

(—a//?) 1 ^ 2 , H c = (Aita 2 /(3) 1 ^ 2 is the thermodynamic 
critical magnetic field. Radially symmetric vortex solu¬ 
tion can be found by solving the following equations 


In the present work, we revisit the problem of the ther¬ 
modynamic stability of the vortex in a superconductor. 
We show that, there is a simple and rigorous relation 
between the magnetic field at the center of a vortex and 
the thermodynamic stability of vortex. Concretely, when 
the magnetic field at the center of vortex is larger than 
the thermodynamic critical field of the material, vortex 
is unstable and the superconductor is of type-1. On the 
contrary, when the magnetic field at the vortex core cen¬ 
ter is smaller than the thermodynamic field, vortex is 
stable and the superconductor is of type-2. The critical 
value, while the field at the core center equals to the ther¬ 
modynamic critical field, represents a boundary between 
two cases. Since the magnetic response of a supercon¬ 
ductor is completely determined by the the stability of 
vortex, our finding can serve as a new criterion for the 
classification of superconductors. 

Vortex solution can always be constructed in the phe¬ 
nomenological Ginzburg-Landau model for superconduc- 
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where k = A/£ is the GL parameter, £ = h/\j2m* |a| 
is the coherence length. Far from the vortex core, the 
order parameter approaches the ground state value and 
the magnetic field vanishes to guarantee the Meissner 
state, 


|^(°o)| = 1, B(oo) = 0. (3) 

The fact that the vortex carries a nontrivial topolog¬ 
ical charge leads to the important conclusion that the 
flux carried by the vortex is quantized in unit of hc/e*. 
Here we consider an isolated vortex line enclosing a sin¬ 
gle flux quantum, which is expected to have the lowest 
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FIG. 1: (Color online) Magnetic field at the center of a vor¬ 
tex 5(0) as function of the Ginzburg-Landau parameter re. 
5(0) = 1 corresponds to re = l/%/2. 


FIG. 2: (Color online) The Gibbs energy difference AG as 
function of the axial magnetic field 5(0). AG = 0 represents 
the boundary between type-l/type-2 superconductors, and 
corresponds to 5(0) = 1. AG —> 0— as 5(0) <C 1 (see text). 


free energy. In this case, asymptotic behaviors of \ip\ and 
B near the center of the vortex can be deduced from the 
quantization of flux f Bpdp = \/2 /k and Eqs. ©, 

\i>\ = a P + (4) 
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where a is the slope of \i/j\, 5(0) is the magnetic field at 
the center of the vortex. With boundary conditions of 
\ij)\ and 5 ([3]) and (©, Eqs. © can be solved. Then 
the vortex solution and a, 5(0) are determined^. Our 
results are more precise and comprehensive than those 
of previous wor k 14 ! 15 . In addition, precise solutions of 
GL vortex can also be obtained by using the iteration 
method developed by Brandt^. 

We note that the difference of vortex solutions origi¬ 
nates from the difference of the GL parameter re. This 
means that 5(0), the magnetic field at the center of vor¬ 
tex, should be a function of the GL parameter re, 

5 (0) = 5(0)(re). (5) 

Further, considering the fact that the flux quantization 
is independent of the GL parameter re, we then conclude 
that the magnetic field at the vortex center, 5 (0), must 
be a monotonic function of the GL parameter re. 

In figure 1, magnetic field at the center of a vortex as 
a monotonic decreasing function of the GL parameter re 
is shown. The critical value re c = l/y/2 corresponds to 
5(0) = 1. The fact that there exists a monotonic rela¬ 
tion between the GL parameter and the magnetic field 
at the vortex core center makes it possible to investigate 


the physical nature of vortex via 5(0) instead of the GL 
parameter re, as we show below. 

In order to study the stability of the vortex, we con¬ 
sider the Gibbs energy difference between the vortex state 
under the thermodynamic critical held and the Meissner 
state AG — G vor t e x(Hc) G Meissner• The same ap¬ 
proach has been used to investigate the stability of vortex 
in a two-component superconductor^. With the defini¬ 
tion of the Gibbs energy density g = / — H • B/47T and 
equation ©> AG can be written in the following form: 
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Note that the Gibbs energy of the vortex state is de¬ 
creasing with the increasing held. If AG > 0, the ma¬ 
terial always behaves as type-1 superconductor since the 
Meissner state plays a dominate role. Once the magnetic 
flux penetrates into the superconductor in the form of 
vortex under certain value of the external held which is 
smaller than the thermodynamic critical held H e < H Cl 
AG < 0 is expected. AG = 0 represents a boundary 
between the stable/unstable vortex state. 

In figure 2, the variation of AG as a function of the 
magnetic held at the vortex center is presented. This 
is a remarkable result. It relates the local physical 
quantity, i.e., magnetic held at the center of a vortex, 
with the global thermodynamic stability of the vortex 
state. When 5(0) > 1, vortex is thermodynamically un¬ 
stable and AG > 0, the superconductor is of type-1; 
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when B(0) < 1, vortex is thermodynamically stable and 
A G 0, the superconductor is of type 2. Zero of AG 
corresponds to -B(O) = 1. 

Asymptotic behavior of AG at B( 0) -C 1 can be esti¬ 
mated as follows. B{ 0) is a monotonic decreasing func¬ 
tion of the GL parameter k and B(0) -C 1 corresponds to 
k 1. In this case, contributions to AG from the vor¬ 
tex core (p < 1 /k) and far range (p > 1) can be ignored. 
In the region 1/k <C p « 1, The approximate solution 
of the magnetic field is B(p) ss V2/«ln(l/p), |^| a; l. 
After some algebra, we find that, as B( 0) C 1 (k> 1), 
AG oc -1 /(v / 2k) -t 0. 

We then derived the relation between the magnetic 
field at the vortex core center and the thermodynamic 
stability of the vortex. Due to the monotonic relation¬ 
ship between the core field B( 0) and the GL parameter n, 
one can use B( 0) instead of n to investigate the stability 
of the vortex. The critical case _B(0) = 1 separate type- 
1 superconductor in which there exists no stable vortex 
(13(0) > 1) from type-2 material, in which there exist 
stable vortex (i3(0) < 1). 

Although the results of present work are based on the 
standard GL model calculations, we believe that, the con¬ 
clusions have general meaning. The configurations of the 
vortices in a charged matter field always bear a certain 
resemblance to each other, including a normal vortex core 
in which the order parameter is suppressed to zero on the 
axis of the core, the restoration of the order parameter 
and the exclusion of magnetic field far from the vortex 
core. Taking into account the fact that the flux carried 
by a vortex is quantized, the difference between vortices 
originates from the difference of the magnetic field at the 
vortex core center. And it is natural that the physical 
nature of vortex is closely related to the quantity 13(0), 
the magnetic field at the center of a vortex. 


The equivalence of our results and conventional meth¬ 
ods which were used to classify superconductors can be 
demonstrated as follows. Physical interpretation of the 
zero of the surface energy is that the thermodynamic field 
is just the critical field at which vortex entry into the 
superconductor. The same case can be extended to the 
vortex geometry, where we set the external field equals to 
the thermodynamic critical field. Besides the emergence 
of a nontrivial topological charge due to the variation of 
the phase of the order parameter around the vortex core, 
when the external field equals to the field at the vortex 
center, i3(0) = 1, AG = 0 can be verified. Moreover, 
we found that, within the capabilities of our numerical 
simulation, no matter what the number of flux quanta is, 
field at the vortex core center 13(0) = 1 and AG = 0 at 
n = 1/ \[2. This result implies that vortices do not inter¬ 
act in this regime. All these results proved qualitatively 
the equivalence of different types of schemes. 

In summary, based on the thermodynamic argument, 
we have shown that the stability of vortex in a super¬ 
conductor, or equally, the classification of the supercon¬ 
ductors according to their magnetic response, has a close 
relation with the magnitude of the magnetic field at the 
vortex center. This finding provides a new criterion for 
the classification of the superconductors. Finally, we note 
that there exist models in high energy physics which are 
mathematically similar to the GL theory for supercon¬ 
ductivity (e.g., the Abelian Higgs modei^ 20 '). Our re¬ 
sults can also be generalized to the investigation of the 
stability of topological vortex in these models. 
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